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We propose, within the effective one body (EOB) approach, a new, resummed, analytical represen- 
tation of the gravitational wave energy flux absorbed by a system of two circularized (nonspinning) 
black holes. This expression is such to be well-behaved in the strong-field, fast motion regime, no- 
tably up to the EOB-defined last unstable orbit. Building conceptually upon the procedure adopted 
to resum the multipolar asymptotic energy flux, we introduce a multiplicative decomposition of the 
multipolar absorbed flux made by three factors: (i) the leading-order contribution, (ii) an "effective 
source" and (iii) a new residual amplitude correction {f>f m ) 21 . In the test-mass limit, we use a 
frequency-domain perturbative approach to accurately compute numerically the horizon-absorbed 
fluxes along a sequence of stable and unstable circular orbits and we extract from them the functions 
pf m . These quantities are then fitted via rational functions. The resulting analytically represented 
test-mass knowledge is then suitably hybridized with lower-order analytical information that is valid 
for any mass ratio. This yields a resummed representation of the absorbed flux for a generic, cir- 
cularized, nonspinning black-hole binary. Our result adds new information to the state-of-the-art 
calculation of the absorbed flux at fractional 5 post-Newtonian order [S. Taylor and E. Poisson, Phys. 
Rev. D 78 084016 (2008)], that is recovered in the weak-field limit approximation by construction. 

PACS numbers: 04.30.Db, 95.30.Sf, 



I. INTRODUCTION 

A ground-based network of interferometric gravita- 
tional wave (GW) detectors is currently being upgraded 
to an improved sensitivity and is expected to detect, 
within a few years, the GW signal emitted during the 
inspiral and merger of compact binaries. The effective 
one body (EOB) approach to the general relativistic two- 
body dynamics [1-12] is the most promising analytical 
tool currently available to model the dynamics and ra- 
diation emitted by coalescing compact binaries, notably 
binary black holes (BBH). The EOB analytical descrip- 
tion of the binary dynamics mainly relies on two build- 
ing blocks: the resummed EOB Hamiltonian H rca \, that 
takes into account conservative effects; and the mechan- 
ical angular momentum loss J- v , also given in resummed 
form, that takes into account nonconservative effects due 
to GW emission (radiation reaction). The current ana- 
lytical construction of both i? rea i and !F V has reached a 
high degree of sophistication for both spinning and non- 
spinning binaries. Several recent comparisons with nu- 
merical relativity simulations [13-17] of inspiralling and 
coalescing relativistic binaries, and with gravitational self 
force calculations [18] have vindicated the EOB approach 
both as a mean to deeply understand the physical out- 
come of such simulations and as a practical tool to build 
accurate analytical template waveforms for data analysis 
purposes [19]. 

In the typical EOB implementation for circularized 
black hole binaries with comparable masses, the radia- 
tion reaction force T v is computed only from the asymp- 



totic GW energy flux. This choice is typically rather 
accurate during the early inspiral phase. On the other 
hand, the cumulative effect during the late inspiral and 
the plunge — typically after the crossing of the last stable 
orbit (LSO) — of horizon absorption can have a nonnegli- 
gible, though small, influence on the phasing. This might 
have some relevance (especially in the spinning case) in 
the ongoing effort in interfacing numerical simulations 
with analytical, EOB-based, predictions. 

The analytical results about horizon absorption of 
GWs are collected in a few papers: in a pioneering study 
in the test-mass limit, Poisson and Sasaki [20] computed 
analytically the leading-order contribution to black-hole 
absorption for a particle on circular orbits around a 
Schwarzschild black hole, providing closed analytical for- 
mulas valid for all multipolcs. This work was then gen- 
eralized by Tagoshi, Mano and Takasugi [21], who con- 
sidered black-hole absorption for a particle in circular 
equatorial orbits around a Kerr black hole. Alvi [22] pro- 
vided the leading order calculation of the growth of mass 
and angular momentum due to horizon absorption for a 
general, spinning binary. 

More recently, Poisson and coworkers have been very 
active in investigating the dynamics of tidally deformed 
black holes [23-26], with several computations (in the 
small- hole/slow- motion approximation) of the tidal heat- 
ing and torquing of a black hole induced by an exter- 
nal tidal environment. In their language, the expres- 
sion "tidal heating" of a black hole refers to the change 
of its mass that occurs as a result of the tidal interac- 
tion; "tidal torquing" refers to the change of its angular 
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momentum. The state-of-the-art analytical information 
that we shall exploit here (in the nonspinning case only) 
is the tidal heating of a black hole of mass m,A created 
by an external body of mass raj as computed by Taylor 
and Poisson [24] for a circularized, spinning, binary in the 
weak-field approximation. The result of Rcf. [24] is ob- 
tained at a higher order of accuracy than what is achieved 
in Ref. [22]. The crucial physical information that can 
be extracted from these analytical studies can be sum- 
marized in two points: absorption effects are fractionally 
more important (i) when the black holes are spinning and 
(ii) when the mass-ratio is large. 

In the comparable mass case, and building on a sugges- 
tion of Price and Whelan [27], the influence of absorption 
effects on the phasing was also investigated by Alvi [22] 
for the first time, using the weak-field approximation and 
at leading order. The conclusion of this study is that 
horizon absorption would change the phasing by no more 
than a tenth of a cycle by the time the binary merge. No- 
tably, this leading-order estimate is much larger than the 
uncertainties of state-of-the-art numerical simulations of 
coalescing black hole binaries [28, 29], suggesting that 
it can be an effect extractable from the numerical data. 
Consistently with this understanding, a recent NR simu- 
lation of a nearly-extremal spinning black hole binary [30] 
has shown a remarkable agreement, during the early in- 
spiral, between the rate of growth of the irreducible mass 
and angular momentum of the black hole horizons com- 
puted numerically and Alvi's analytical prediction. This 
numerical result suggests that it is recommendable to in- 
clude horizon absorption effects in the EOB-based ana- 
lytical modeling of coalescing black hole binaries. To our 
knowledge, absorption terms have been included in the 
EOB model (though in nonresummed form) within the 
context of an extreme-mass-ratio (spinning) inspiralling 
binary, a study that aimed at determining the suitability 
of EOB-based template waveforms for LISA data analy- 
sis purposes [31, 32]. 

For noncircularizcd binaries the effect of black-hole 
absorption is more relevant [33-35]. For example, in 
the test-mass limit (on Schwarzschild) absorption effects 
were thoroughly investigated by Martel [34] (and more 
recently in Ref. [36]): he considered both eccentric and 
parabolic orbits and found that the absorbed flux is non- 
ncgliblc (larger than 1% of the asymptotic flux) for or- 
bits whose periastron is smaller than 5M. On the ba- 
sis of these test-mass results, one expects that horizon 
fluxes should have some relevance (even in the nonspin- 
ning case) also for comparable mass binaries with an 
analogous dynamical setup, for example when a zoom- 
whirl behavior is present [37-40]. 

The main purpose of this paper is to provide an analyt- 
ical expression of the (multipolar) horizon-absorbed flux 
for circularized, nonspinning, binaries within the EOB 
framework that is valid also in the strong-field, high- 
velocity regime, i.e. possibly up to the EOB-defincd last 
unstable orbit. This result builds upon three different 
pieces of knowledge: (i) the leading-order and next-to- 



lcading order analytical results for the multipolar horizon 
flux of Refs. [20, 24]; (ii) the accurate numerical computa- 
tion of the horizon fluxes for a test-particle orbiting along 
stable and unstable circular orbits of a Schwarzschild 
black-hole; (iii) the idea of suitably factorizing the an- 
alytical representation of the horizon flux, analogously 
to what was proposed (actually, at the waveform level) 
in Rcfs. [5, 41]. In addition, the numerical computation 
of the horizon-absorbed flux for unstable circular orbits 
of Schwarzschild is presented here for the first time. 

The paper is organized as follows: Sec. II reviews the 
next-to- leading order analytical results of Refs. [20, 24]. 
In Sec. Ill we first present the structure of the factorized 
resummation of the asymptotic multipolar energy flux 
for circularized binaries within the EOB approach and 
then wc introduce our, conceptually analogous, new, rc- 
summcd representation of the multipolar horizon flux. In 
Sec. IV we calculate the residual amplitude corrections 
Ptm fr° m the numerically computed horizon flux in the 
test-mass limit and in Sec. V we show how to use them 
together with the analytical results of Sec. Ill so as to ob- 
tain an explicit analytical expression of the horizon flux 
that holds for any mass ratio. We collect our conclusions 
in Sec. VI. We use geomctrized units with G = c = 1. 

II. ABSORBED ENERGY FLUX: 
PN-EXPANDED ANALYTICAL RESULTS 

We start by reviewing the analytical results of 
Rcfs. [20, 24]. We consider a system made of two nonspin- 
ning black holes of masses vtla and tub (with total mass 
M = rriA + tub) in circular orbit. In the following we 
shall use the symmetric mass ratio parameter v = /x/M, 
where /i = m^ws/M. This parameter varies between 
v = (test- mass) and v = 1/4 (equal- mass) . Here, the 
"test- mass" limit is defined by the condition ms <C rriA, 
i.e. M ~ uia becomes the mass of the larger black-hole, 
/i ~ rris that of the smaller one and v ~ mB/'niA- The 
total energy flux emitted by the system can be written 
as a multipolar expansion of the form 

1=2 m=l 

where each multipolar contribution 1 F% m is given by the 
sum of the flux emitted at infinity, and the "horizon" 
absorbed flux, Ff m , 

Ftm = i^m + Ff m . (2) 

In the test-particle limit the absorbed flux is the energy 
flux through the event horizon of the large mass. In the 



1 Our definition assumes that the +m and —m modes are summed 
together. 
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case of two Schwarzschild black holes, it is the sum of the 
fluxes through the two horizons 2 . 

For a test-mass on Schwarzschild circular orbits the 
functions have been recently computed by Fujita [43] 
up to the 14 post-Newtonian (PN) order. This achieve- 
ment builds upon a systematic method that treats black- 
hole perturbations [44] and it improves on previous state- 
of-the-art 5.5PN-accurate calculation [45]. For the gen- 
eral ^-dependent case, the PN knowledge of F^ n has been 
pushed up to 3.5PN order [46-49]. 

As mentioned in the introduction, much less analytical 
information is available for Ff m . The leading order (LO) 
contributions (for each multipole) in the test-particle 
limit (with the mass-ratio fi/M = v <C 1) were computed 
long ago in closed form by Poisson and Sasaki [20] . Com- 
bining their Eqs. (5.7) and (5.13)-(5.18), one can rewrite 
each multipolar contribution in compact form as 
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where x = (A/.!?) 2 / 3 is the PN frequency parameter (J? 
is the system orbital frequency) and e denotes the parity 
of the flux, i.e., even (e = 0) for mass-generated mul- 
tipoles and odd (e = 1) for current-generated ones. In 
the circular case, e is equal to the parity of the sum 
I + to : e = tt(£ + to), that is e = when I + to is even 
and e = 1 when £ + to is odd. The quantities oYg m and 
_iYf m are the spherical harmonics of spin- weight s = 
and s = — 1, that are given explicitly in Eqs. (5.7) and 
(5.13) of Ref. [20], while 
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When the masses are comparable (and possibly equal), 
the state-of-the-art analytical knowledge is due to Tay- 
lor and Poisson [24]. As in the test-mass limit, the cal- 
culation is based on the weak- field approximation but, 
contrary to the result of Ref. [20] , it is complete at next- 
to-leading order (NLO), i.e. 1PN fractional accuracy. 

Let us now review the results Ref. [24] and recast them 
in a form that is useful for our purpose. Focusing on the 
object labeled by A, the NLO (i.e., 10PN accurate) "rate 
at which the black hole acquires mass by tidal heating" 
is expressed as the sum of a £ = to = 2 and a £ = 2, 



2 Note that our definition of F? in the test-mass limit is neglect- 
ing the flux absorbed by the black hole of smaller mass that is 
in any case smaller by the fourth-power of the mass ratio. It 
is however remarkable that the quadrupolar contribution to this 
quantity has been explicitly computed by Poisson in closed form, 
see Eq. (8.52) of Ref. [42]. If the need arises, it can just be added 
to Eq. (2). 



to = 1 contributions, which explicitly read (see Eq. (9.1) 
of [24]) 

•22 16 6 ? j? ab 
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where f is the orbital separation in harmonic coordinates. 
The quantities £ a b and B a b are the electric and magnetic 
tidal moments in the black hole frame and the overdot 
refers to time derivative in the same frame; we used the 
relations mentioned after Eqs. (9.2)-(9.3) of [24]) to write 
explicitly the right-hand sides. Summing together these 
two pieces gives the final, fractionally lPN-accurate, con- 
tribution to the expression for tidal heating given by 
Eq. (9.4) of Ref. [24]. The total multipolar horizon fluxes 
of the binary system are then given by 



Fg{f- v)=m A 2 +m% 2 , 
Fg(f- v) = m\ x 



which become 
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where we made now explicit the dependence on the sym- 
metric mass ratio parameter v. One then wants to 
express the harmonic radius r in terms of the gauge- 
invariant frequency parameter x = (MJ7) 2 / 3 . At the 
needed 1PN accuracy, this relation reads [50] 



M 

f 



1 



1 



) x + 0(x 2 ) 



(12) 



so to get the gauge-invariant results 
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Note that the leading-order contributions 
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reduce to Eq. (3) (for I = 2, m = 1,2) in the v 
limit. Equations (13)-(14) are basically the only ana- 
lytical information present in the literature about hori- 
zon absorption for nonspinning black-hole binaries in the 
comparable mass case 3 . 



III. RESUMMED REPRESENTATION OF THE 
ABSORBED ENERGY FLUX 

The purpose of this section is to obtain a physically 
motivated analytical expression of the absorbed energy 
flux for circularized, black-hole binaries (for any mass 
ratio v) that has the property to remain valid also away 
from the weak-field, slow-motion assumption that is at 
the basis of the derivation of Eqs. (13)-(14). More pre- 
cisely, we would like to obtain, by a physically motivated 
guess within the EOB framework, a resummed analytical 
expression that is well behaved in the strong-field, fast- 
motion regime, possibly up to the last unstable circular 
orbit. To to so, we will introduce a factorized, analytical 
model of the absorbed multipolar flux that is inspired by 
the procedure followed in Ref. [5] to treat the analogous 
problem in the case of asymptotic fluxes. 



A. Factorization of the multipolar asymptotic flux 

Let us first recall the prescription that is used to com- 
pute the resummed asymptotic energy flux (in a factor- 
ized form) for circularized binaries in the EOB approach. 
For circular orbits, the dynamics is typically parameter- 
ized in terms of the gauge-invariant frequency parameter 
x = (MQ) 2 / 3 introduced above. The conservative dy- 
namics is described by the real EOB Hamiltonian 



H iesl = MJl + 2v Lffeff-1 



(17) 



where H c g = H c s/fi is the (reduced) effective Hamilto- 
nian. When v — s- 0, H c g reduces to the usual conserved 
energy of a test-mass fx in a Schwarzschild background 
of mass M. The explicit expression of H g as a function 
of the frequency parameter x cannot be written in closed 
form, but rather in parameterized form in terms of the 
EOB inverse radius parameter u = M/r, that is as 



ffeff = V 'Mu)(l + J 2 U 2 ) 



(18) 



where A(u) is the EOB radial potential and j is the 
dimensionlcss angular momentum along circular orbits. 



3 Notably, in Ref. [23] Poisson was able to compute also the oc- 
tupole contributions to the metric of a distorted black hole. This 
means that in the literature there is already enough information 
to compute at leading order the I = 3 black-hole absorption part 
by following the same procedure of Ref. [24] . 



The PN-expansion of A(u) at 3PN accuracy is 

'94 41 
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(19) 

but then we use it in a Pade resummed form, that allows 
for the presence of an "effective" horizon, an effective 
"light-ring" (i.e., corresponding to the last unstable orbit) 
and of the last stable orbit (LSO). Since we work at 3PN 
accuracy, we use a (1,3) Pade representation 



A{u) = Pi [A Taylor ] = 



no + n\u 



(20) 



da + d\u + d2U 2 + dzu 3 

The circular orbits in the EOB formalism are determined 
by the condition d u {A(u)[l + j 2 u 2 ]} = which leads to 
the following parametric representation of the squared 
angular momentum 



f(u) = - 



A'(u) 



(21) 



(u 2 A(u))'' 

where the prime stands for d/du. Inserting this u- 
parametric representation of j 2 in Eq. (18) defines the 
it-parametric representation of the effective Hamiltonian 
H c g(u). One can then obtain H e g as a function of x 
by eliminating u between H e g(u) and the corresponding 
it-parametric representation of the frequency parameter 
x = (ML?) 2 / 3 obtained by the angular Hamilton equation 
of motion in the circular case 

M n(u) = ^^= MA ^ u)u2 . (22) 
M 9j H ieal H eS 

The use of the variable x to parameterize the adiabatic 
circular dynamics is the standard choice when one focuses 
on stable circular orbits [3, 5] (i.e., orbits corresponding 
to the local minimum of the effective potential A(u)(l + 
j 2 u 2 ) [1, 51] for a given value of j). In this case, the 
factorization of the waveform introduced in [5] allows us 
to factorize each multipole of the asymptotic energy flux 
as 
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where F^' e \x; v) is the Newtonian contribution for 

circular orbits, (x; v) symbolically represents all 

higher-order relativistic corrections and e is the par- 
ity of the sum £ + m. To be precise, the quantities 
F^ (x; v) cx x 3+i+e are the fluxes computed from 
the Newtonian contribution to the circular gravitational 
waveform, as given by Eq. (4) of Ref. [5]. The rel- 

ativistic correction factor is written in resummed form 
as 

2 



F, 



(00,£) 



where S^g is the effective source of the field, 

S^(x) = Hctt(x) £ + m even, 

(x) = j(x) = x 1/2 j(x) £ + m odd. 



(24) 
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The quantity T( m is the "tail factor" that takes into ac- 
count the infinite number of "leading logarithms" enter- 
ing the transfer function between the near zone multipo- 
lar wave and the far-zone one due to tail effects. This 
factor is written as [5] 



1 - 2ik) 



T(l + 1) 



e 7r/c e 2ifelog(2fer ) 
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where k = m(H rcll i/M)x 3 ^ 2 , k = (m/M)x 3 ^ 2 and r = 
2/y/e [52]. Finally, pf^ are the residual amplitude correc- 
tions that are essentially given as Taylor series (modulo 
log{x) n terms) of the form p^ n {x; v) = l+x+x 2 +x 3 + . . . 
with ^-dependent coefficients. In the test-mass limit, 
the PimS> are currently analytically known at 14PN or- 
der [43]. In the general, ^-dependent case the available 
analytical information is more limited, nonetheless it gets 
up to 3PN order [5]. This independent information is 
then "hybridized" with the test-mass limit results (at 
5PN-accuracy only) to compute the so-called 3 +2 PN ap- 
proximation that has been used recently in several com- 
parisons between EOB predictions and numerical simula- 
tions of nonspinning inspiralling and coalescing compact 
binaries [15-17, 41, 53]. Moreover, in the test-mass limit 
the p^j's have also been computed numerically, using 
black-hole perturbation theory, for stable [5] and unsta- 
ble [54] circular orbits. 

When one wants to consider also the unstable circu- 
lar orbits predicted by the EOB conservative dynam- 
ics [1, 51] (i.e., orbits corresponding to the local maxi- 
mum of the effective potential A(u)(l +j 2 u 2 ) for a given 
value of j), as we are going to do here, a complication 
occurs, since the variable x is found to be ill-behaved for 
these orbits. When expressed as a function of u, x does 
not grow monotonically (as in the test-mass limit, where 
it is simply x = u), but it has a maximum and then 
decreases, to go to zero at the EOB last unstable orbit 
(or EOB-defined "light-ring", that in the adiabatic ap- 
proximation is defined as the maximum of the function 

u 2 A{u)). The effect of this on the fluxes F^' e \x; v) 
is that they grow monotonically until the maximum of 
x and then turn back (i.e., they are not functions any- 
more). This behavior is obviously qualitatively incom- 
patible with the test-mass result, where the asymptotic 
fluxes are always growing monotonically until they di- 
verge at the light-ring position, x = 1/3. 

To meaningfully compute a resummed representation 
of the energy flux along unstable orbits we need first to 
use a different, well-defined, frequency parameter. A sim- 
ple possibility is given by the following variable 



(if real i2) 
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M 



2/3 
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which does not vanish at the EOB light-ring as a func- 
tion of u and it is continuously connected to the standard 



frequency parameter x both in the test-mass and New- 
tonian limits. The notation x that we use here reminds 
us of the variable k in Eq. (27), that can be written as 



mx 



3/2 



More importantly, the use of x suggests a 
slightly different factorization of the resummed energy 
flux. To derive it, we take Eq. (23) and we first make it 
explicit the dependence of x on x, i.e. 
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where x[x) symbolically represents the inverse of 
Eq. (28). Note that in the tail factor the dependence 
on x is left explicit because of Eq. (27). We then replace 
the function x(x) in F^^ (x(x); isj and p'^ n (x(x); v) 
by 



Z3PN(£) = Xf(x) 
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that is obtained by inverting the 3PN-accurate expan- 
sion of Eq. (28), and, as a last step, we crucially factor- 
ize out from the resulting modified Eq. (29) the Newto- 
nian prefactor considered now as an explicit function of 
x, F^' e \x) oc x e+3+e . It follows that a new residual am- 
plitude corrections /5^(i; v) is defined as the following 
expansion in powers of x 



Pirn 



(x; v) = Taylor 
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{/( 



(31) 

where the notation Taylor^ {...} indicates the Taylor ex- 
pansion of the argument inside the curly brackets at order 
N. Here N is the highest power present in the Taylor- 
expanded p^ix; v) in Eq. (24). At the end, a new, dif- 
ferent, ^-dependent, expression for the asymptotic mul- 
tipolar flux that is well-behaved along both stable and 
unstable circular orbits is 
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Note that the z^-dependent coefficients of the Taylor- 
expanded function p^ n (x; v) are different from the cor- 
responding ones of the p < ^ n {x] v). On the contrary, in 
the test-mass limit (where x — x) one has p^{x; 0) = 



pf m (x- 0) (and thus F^\x; 0) = F^>(x; 0)). Since 
we shall need it below, we explicitly write the expression 
for p^(x; v) expanded formally to 3PN (in the variable 
x) accuracy 
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where eulerlog m (x) = 7 + log 2 + — logx + logm. Evi- 
dently Eqs. (31)-(33) define only one among the many 
possible Taylor-expanded representations of p^(x] v). 
For example, to quote one alternative possibility, we 
could have used the Taylor-expansion of Eq. (28) to ex- 
press x as a function of x and then use again the non- 
expanded implicit relation x(x) to represent the result 
as a function of x. Whatever one does for representing 
the relativistic correction, the crucial new element that 
makes the flux well behaved also below the LSO is the 
factorization of the Newtonian contribution as an explicit 
function of x and not of x. Note that, for simplicity, 
we did not include in Eq. (33) also the 4PN and 5PN 
test-mass terms that can be used to compute the hybrid 
3 +2 PN approximation to p^(£; v). They can be found 
in Eq. (50) of Rcf. [5] and added to Eq. (33) as i 4 and 
x 5 terms. 



B. Factorization of the multipolar horizon flux 

We have seen in the previous section that after factor- 
ization all information about the asymptotic GW ampli- 
tude and asymptotic GW energy flux is contained in the 
p'fan (or pfj^) functions that are therefore the only real 
unknowns of the problem (for a given factorization). Let 
us focus now on the horizon-absorbed flux. The strat- 
egy that we shall follow here is to introduce a resummcd 
factorization procedure for the multipolar horizon fluxes 
that is analogous to the one used above for the asymp- 
totic fluxes. The main idea is to factorize out from the 
PN-expanded representation of the horizon fluxes, given 
by Eqs. (3), (13) and (14), both the leading-order con- 
tributions (analogous to the Newtonian prefactor in the 
asymptotic fluxes) and the effective source so to eventu- 
ally introduce some residual (horizon) amplitude correc- 
tions pf m (x; v) of the form 1+.t+.t 2 + . . . (or correspond- 
ingly pf m (x; v) = 1 + x + x 2 + . . . ) that take into account 
the remaining PN corrections. As in the previous sec- 
tion, we first present a factorization that is well-behaved 
for stable circular orbits only. We then start from this 
result to obtain a different factorization of the horizon 
flux where the leading-order prefactor is an explicit func- 
tion of x and so it is meaningfully defined for stable and 
unstable circular orbits. 

We write each multipolar contribution to the horizon 



flux as 

F ( " n e \x; = F^°' e \x; v)F^\x; „), (34) 

where F^" 3 {x; v) indicates the leading order contri- 
bution to the absorbed fluxes (given by Eqs. (15)-(16) 
for £ = 2 and Eq. (3) for the other multipoles 4 ) and 

Fg m ' '(x; v) accounts for all higher-order relativistic cor- 

rections. We then factorize F^ ' '(x; v) as 

e) (*; ") = (Scff(z)) 2 (pL(*\ ")f ■ (35) 

The physical reason for factorizing the effective source 
in Ff m ' {x\ v) is the same one done in Rcf. [5] for 

Ftm \ x 'i v )- The idea in Ref. [5] was motivated by the 
form of the equation satisfied by each partial wave in 
the circular test-mass limit with a source made up of the 
important dynamical constituents of (test-particle) dy- 
namics, i.e. its energy and angular momentum (see dis- 
cussion below Eq. (8) of [5]). Note that, as it is the case 
with F^'^ (x; v), the factorization of the source for each 
multipole actually means factorizing out the singular be- 
havior of the absorbed energy flux at the (EOB-dcfined) 
light-ring 5 . Notably, Ref. [54] has shown that, in the 
test-mass limit, the numerical p^Ax; 0) functions appear 
to be regular when x is close to 1/3. This a posteriori 
confirms the intuition that the singular behavior of the 
flux at the light-ring is only due to the squared source. 
Similarly, in proposing here the factorization (35), we are 
assuming the same singular behavior also for the horizon 
fluxes, so that we expect the functions pf m (x; 0) to be 
regular when x is close to 1/3. In Sec. IV we will see 
that this intuition is supported by the behavior of the 
"exact" test- mass pf m {x; 0) computed numerically 6 . 



4 Note that in this context v is always assumed to be the symmetric 
mass ratio v = m^ms/M 2 , even if Eq. (3) is obtained in the 
test-particle limit. In doing so we are in fact neglecting the (yet 
unknown) polynomials in v, of the form 1 + v+u 2 + . . . , that have 
to appear in the LO prcfactors when the ^-dependent calculation 
of Eqs. (15)-(16) is extended to the other multipoles. 

5 This idea is similar to the suggestion of Ref. [55] of factoring out 
a pole from the energy flux. 

6 In this respect, note that the absorbed flux in the small-hole ap- 
proximation computed by Poisson in Ref. [42], Eq. (8.52), has a 
rather different structure since the singular behavior is propor- 
tional to the fourth power of the source. 
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Let us now compute the explicit analytical expressions 
of the pf m (x; v)'s. Evidently, since (1PN) fractional cor- 
rections have been analytically computed only for the 
t = ui = 2 PN-expanded horizon flux, Eq. (13), one can 
obtain them only for p^(a;; v) and only at 1PN fractional 
accuracy, while pf m {x; v) = 1 (i.e., OPN fractional accu- 
racy) for all other multipoles. By factoring out of Eq. (13) 
the effective source, H e g(x; v), and the leading-order 
contribution, F 2 f LO (x; v) = (32/5)i/ 2 (l - Av + 2v 2 )x 9 , 
taking the fourth root and then expanding this latter at 
1PN accuracy, we obtain 



4-2lv + 27v 2 



8v c 



4(l-4v + 2v 2 ) 



x + 0{x 2 ). (36) 



The combination of Eq. (35) and Eq. (36) gives an ana- 
lytical expression of the multipolar horizon flux that, as 
it was the case for Eq. (23) is ill-behaved along unstable 
orbits. We need then to apply the same procedure that 
lead to Eq. (32) for the asymptotic flux. We first make 
explicit the dependence on x in the horizon flux 



im 



(x(x) 



Mi 



(x(x); v) 
(S^(x(x))) (pf m (x(x); u) 



(37) 

then we replace in both Fp^° (x(x); v) and 
pf m (x(x); v) the function x(x) with X3pn(x), Eq. (30), 

'(x) as an 



(Hlo, 



and in the end we factorize in front F t 
explicit function of x. Since this factor is proportional to 
x 2f+5+e , this operation yields a new residual amplitude 
correction as the following expansion in powers of x 



Ptm& v ) = Taylor^ |/(x) 



( 2 , + 5 +£ )/( 2 ,) pfm[x3pNl 



i(*)]} > 
^ (38) 

where N is now the highest power present in the Taylor- 
expanded pf m {x; v) of Eq. (35) above. Eventually the 
new, x-dependent, expression for the horizon multipolar 
flux that is well-behaved along stable and unstable orbits 



x (S^(x(x))) 2 [^(i; v) 
In the I = m = 2 case we have 
P22 (x; v) = 1 + 



2/ 



(39) 



4-18^ + 15^ 2 -2^ 3 : 



4(1 - 4v + 2v 2 ) 



x + C(x 2 ), (40) 



while it is pe m (x; v) = 1 for all other multipoles. For 
v = 0, where x = x and thus pf m {x; 0) = pf m {x; 0)), this 
equation reduces to the very simple result /o|^(x; 0) = 
1 + x + 0(x 2 ). Note also that the dependence on v 
is relatively mild, since the coefficient of x in Eq. (40) 
varies between 1, when v = 0, and 13/16 ~ 0.8125, when 



0.18 
0.16 
0.14 - 
0.12 - 



tt( 0.08 - 



-F H /v 2 (4»x = 8) 




FIG. 1. Horizon-absorbed energy flux for a point-particle 
orbiting along stable and unstable circular orbits of the 
Schwarzschild spacetime. The total flux (summed up to 



is contrasted with the quadrupolar £ 



flux. The same quantities are plotted on a linear (top) and 
logarithmic (bottom) vertical scale. Here x — 1/r, where r is 
the Schwarzschild radial coordinate. The light-ring is located 
at x = 1/3. 



v = 1/4. This corresponds to less than a 5% fractional 
difference at the test-mass light-ring position, x = 1/3 
between v = and v — 1/4. In the next Section we shall 
check the consistency of the 1PN analytical prediction 
and its accuracy with respect to the /of?,(x; 0) computed 
numerically 



IV. NUMERICAL COMPUTATION OF 



0). 



Since the highest PN analytical information available 
for the pf m 's is given by the lPN-accurate result of 
Eq. (40), that is limited to the £ = m = 2 multipolc, 
the first question that one wants to address is about its 
accuracy in the strong-field, fast- velocity regime. In the 
test-mass limit this question can be answered exhaus- 
tively by comparing Eq. (40) with the numerical /of?>(x; 0) 
that can be extracted from the numerical computation 
of the GW absorbed flux, in a way similar to what was 
done in Ref. [5] for the asymptotic fluxes. To this aim, 
we computed the horizon absorbed energy flux for a se- 
lected sample of circular orbits belonging to both the 
stable and unstable branch. Although the horizon flux 
has been computed in the past for stable circular or- 
bits [34, 56, 57], we are not aware of any published results 
for unstable orbits. We computed the horizon fluxes us- 
ing the gravitational self-force (GSF) code developed in 
Ref. [58] , which is a frequency domain variant of the orig- 
inal self- force code of Ref. [57] . The GSF code used com- 
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FIG. 2. The numerical pf m (x; 0) computed from the numerical fluxes using Eq. (41). The dashed (red online) lines indicate 
the outcome of the fit using the rational function given in Eq. (44). 



putes the 0(/i/M) metric perturbation. The flux is then 
obtained in terms of surface integrals of Weyl scalars eval- 
uated at the event horizon of the large mass. For circu- 
lar orbits, the expressions for the horizon and asymptotic 
fluxes go back to Teukolsky's work in the 1970s [59, 60]. 
We consider orbits with radii varying from r m i n = 3.1 to 



r max = 30, that are spaced by Ar = 0.1 for r <G [3.1, 10] 
and by Ar = 1 for r <G [10,30]. We have computed the 
numerical F^/v 2 up to £ m &x = 8: Fig. 1 compares the 
total flux (summed up to ^ max = 8) with the £ = m = 2 
contribution. The frequency domain approach allows us 
to obtain the F^ uin / V 2 's very accurately. In particu- 
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TABLE I. Values of the coefficients of the fitting function Eq. (44) for multipoles up to I — 7. 



£ 


J7l 


ni m 

1 "1 


' "2 




n| m 




sr 


2 


2 


- , 7801 6 


7 022Q1 


7096 . 


-33 06646 


-4.71422 




2 


1 


-1 89997 


49519 


5 98688 


-17 15286 


-2.48118 


92669 


3 


3 


-2 36556 


1.73143 


34 37095 


-124 89266 


-3 50205 


1.87045 


3 


2 


-1.47910 


-0.37517 


20.54194 


-56 00960 


-2.31621 


16676 


3 


1 


-1 50854 


0.24240 


4.47148 


-9 36740 


-3.11918 


2.29453 


4 


4 


-1.96477 


1.77455 


41 82985 


-155 50256 


-3.11767 


0.77264 


4 


3 


-0 88873 


-9 35298 

Zj . (J Zj _ O 


38 83989 


-1 1 8 591 1 3 


-1.84937 


-9 03044 


4 


2 


-0 35899 


-0.14473 


1 7 55970 


-98 91 807 

ZjO. -LOU 1 


-1 79350 


-0 0041 3 


4 


1 


-1 86354 


1 6555 


3.23144 


-8.41344 


-2.76941 


1 49953 


5 


5 


-2 01738 


3.91457 


35 23863 


-150 27033 


-3 12766 


87106 


g 


4 


-0 89483 


-1 79998 


44 98785 


-154 35735 

±<J^L. tJ'J I >J'.J 


-1 88958 


-2 54359 


5 


3 


3 96179 


-1 42399 


50 85186 


-23 44662 

Z » J . I 1 V f V ■ z 


2 58708 


-5 91 327 


5 


2 


-0 79955 


-1 18556 


1 3 74389 


-34 0331 3 


-1.76133 


-0.62741 


5 


1 


-1.68170 


0.13047 


2.74892 


-6 05858 


-3.07511 


2.19172 


g 


6 


-2 27965 


6 83094 


21 35950 


-125 03545 


-3.33171 


1 53708 


6 


5 


-1 35000 


1.11200 


36 38329 


-148 13578 


-2 32261 


-1 50432 


g 


4 


17.75737 


-4 90987 


160 53963 


101 19617 


1 6 39948 


-21 05465 




3 


4 39780 


1.01247 


43 69830 


9 99466 


3 37649 


-9 661 69 


g 


2 


-0 52866 


-0 68194 


1 1 73839 


-24.81141 


-1.85746 


-0.15371 


6 


1 


-1 88052 


1 7089 


2 21690 

Zj ■ Z. 1 V ' i.' V ' 


-5 55880 


-2 83756 


1 68356 


7 


7 


-2.60076 


9.77862 


5.23594 


-92.25665 


-3.59794 


2.37691 


7 


6 


-1.90285 


4.69608 


20.96354 


-120.14340 


-2.83491 


-0.03207 


7 


5 


159.88998 


-25.98598 


1348.46529 


1881.03562 


157.35115 


-161.73362 


7 


4 


19.50198 


10.68544 


143.28131 


252.11345 


18.36505 


-4.37623 


7 


3 


2.53654 


-2.61813 


33.54445 


-50.19582 


1.23298 


-5.23258 


7 


2 


-0.71568 


-1.37276 


9.79925 


-25.74536 


-1.70002 


-0.80200 


7 


1 


-1.76491 


0.15087 


1.99348 


-4.43712 


-3.04998 


2.14649 



lar, the fractional uncertainty is extremely low, i.e. on 
the order of 10 -10 or smaller for all multipoles and for 
strong-field orbits (say r < 10); it is at most of order 
lCT 2 -lCr 3 for the I = 8, m = (2,3) multipoles when 
r = 30. This apparent "inaccuracy" is due to the fact 
that beyond r ~ 15, the magnitude of (£,m) modes of 
horizon flux is comparable to or less than machine ac- 
curacy. So the concept of "fractional accuracy" becomes 
numerically meaningless. For example, the I = 8, m = 1 
multipolar contribution is (~ 10~ 33 ). 

The pf m (x\ 0) are computed from the horizon numer- 
ical fluxes F^ am simply by 



where in the test-mass limit one has 



Pem( x > °) 



F, 



(in 



y4^°' e) (*; o) 



(41) 



1 - 2x 



sfl -3x' 
Vl — 3a; 



(42) 
(43) 



The result of this computation is shown in Fig. 2 for mul- 
tipoles up to I = 7. An important feature highlighted 
by the plot is that the value of pf m increases by a con- 
siderable factor (between 3 and 4) when x varies from 
1/30 to 1/3.1. Such a behavior is rather different from 
the quasi-linear (decreasing) trend of the numerical pf^ 
found in Refs. [5, 54] (see Fig. 3 and 4 of [5] and Fig. 7 
of [54]). This fact shows that the 1PN approximation 
of P22 will be rather rough and more analytical informa- 
tion is needed to reproduce the correct behavior of the 
pf m {x; 0) of Fig. 2 towards the light-ring (see also the 
discussion around Fig. 3 below). 

In the absence of high-order PN results for the horizon 
fluxes that allowed us to push the analytical knowledge 
of pf m to higher order, for the moment we content our- 
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TABLE II. Numerical evaluation of the 1PN coefficient 4 m in 
Eq. (44) from the fit to numerical data. Note that c 22 differs 
from its analytically predicted value (i.e., c 22 = 1 as obtained 
from the v = specification of Eq. (40)), by a 6.6% due to 
the global nature of the fit. 



1 


m 


c{ m 


2 


1 


0.58121 


2 


2 


0.93406 


3 


1 


1.61064 


3 


2 


0.83711 


3 


3 


1.13649 


4 


1 


0.90588 


4 


2 


1.43458 


4 


3 


0.96063 


4 


4 


1.15290 


5 


1 


1.39341 


5 


2 


0.96178 


5 


3 


1.37472 


5 


4 


0.99476 


5 


5 


1.11029 


6 


1 


0.95704 


6 


2 


1.32879 


6 


3 


1.02131 


6 


4 


1.43489 


6 


5 


0.97261 


6 


6 


1.05206 


7 


1 


1.28508 


7 


2 


0.98434 


7 


3 


1.30356 


7 


4 


1.13693 


7 


5 


2.53883 


7 


6 


0.93206 


7 


7 


0.99718 



selves to have an effective analytical representation of the 
numerical results in terms of simple functions. A conve- 
nient procedure to do so is to fit the numerical data with 
a suitable function that tends to 1 when x — > 0. Wc 
found that the following rational function 

p ^ x) = — 1+^+4^ (44) 

works well for all multipolcs (and especially for the £ = 
m = 2 one, which presents a characteristic "bump" close 
to the light ring, similarly to the pi m {x; 0)'s [54]). The 
coefficients are listed in Table I for multipoles up to £ = 7. 
One notes that their values are (on average) relatively 
small due to the choice of a rational function to perform 
the fit. Had we used instead a polynomial expression 
(with the same number of free parameters) we would have 
obtained much larger numerical values. 



From Eq. (40) we have p 22 {x; 0) = 1 + x + 0(x 2 ). It is 
then interesting to compute a PN approximation of the 
fit to , in order to have a concrete idea of the reliability 
of the ansatz (44). By expanding Eq. (44) around x = 
we find 

PtmW = l+c[ m x+c i 2 m x 2 + c e 3 m x 3 + c{ m x 4 +0{x 5 ), (45) 

where cf" = n[ m - d{ m , c e 2 m = {d[ m f - d{ m - d[ m n[ m + 
n 2 and the corresponding dependence of the other cf m 
on nf m and d\ m . The numerical values of c[ m are listed in 
Table II for completeness. First of all, one immediately 
notices that the c[ m are always of order unity. The small- 
est value is cf = 0.58121, and the largest c[ 5 = 2.53883, 
where the fit does not seem to be visually very accurate 
when x — > 0. Inspection of Table II shows that the value 1 
is a good approximation for several (actually, almost all) 
multipolcs. In particular, one has c 22 = 0.93 4 06, which 
is smaller only by a 6.6% with respect to the analytically 
predicted value. This, rather small discrepancy is due to 
the global nature of the fit. Note in fact that this agree- 
ment can be further improved by including also a a; 5 -type 
term in the numerator of Eq. (44), so to allow for more 
flexibility. In this case, the coefficients of the fit take the 
following numerical values {ri\ , ri2, n.3, n^, 715 , d\ , di) = 
(-4.59, 6.01, 26.95, -178.53, 276.68, -5.60, 8.05), which 
yield c 22 = 1.01590. However, the fact that the 7x4 and 
n.5 coefficients become rather large (and relatively close) 
is probably an indication that we are allowing for exces- 
sive freedom in the analytical model. For this reason we 
prefer to work with a 4th-order polynomial at the nu- 
merator even if the results is slightly less accurate when 
compared to the 1PN analytical value. 

An alternative possibility is to impose the analyti- 
cal 1PN value as a constraint in the fit, i.e. by ask- 
ing that d 22 = n 22 — 1. In doing so, the fit im- 
proves (though at a level indistinguishable on the scale 
of Fig. 2, see Fig. 3) and the low-frequency limit is re- 
covered correctly by construction. In this case one gets 
the following numerical values (ni, 712, «3, 714, n^, d%) = 
(-3.51327, 5.28819, 10.16670, -52.90979, 5.01393). Fig- 
ure 3 focuses on the I — m = 2 mode only and in partic- 
ular it illustrates the performance of such a constrained 
fit. The exact data (black solid line) are contrasted 
with the analytical 1PN result (thick, dotted line), to 
their constrained fit (dashed red line) via Eq. (44), and 
to the 2PN-accurate (dashed line) and 4PN-accurate 
(dash-dotted line) approximations of this latter. In this 
case, the 2PN, 3PN and 4PN coefficients are given by 
cf = 4.78752, cf = 26.76014 and cf = 43.86148. As 
mentioned above, the important information conveyed by 
this comparison is that the 1PN analytical result gives 
only a rough approximation to the numerical data, since 
it differs from them by about a factor of 3 when x — > 1/3. 
Interestingly, the plot shows that one needs to expand 
Eq. (44) to (at least) 4PN (effective) order to have a 
good averaged approximation to the full data. This re- 
sult also suggests that one would probably need at least 
a 4PN analytical calculation (beyond the leading order 
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FIG. 3. Test-mass limit: investigating the properties of 
pxiix; 0). The "exact" numerical result is contrasted with: 
(i) the lPN-accurate analytical prediction; (ii) the fit of the 
numerical data by the rational function given by Eq. (44) with 



the lPN-constraint d 1 



nf 



1; (hi) the 2PN-accurate and 



4PN-accurate expansions of this fitting function. 



x 9 contribution) of the absorbed flux (i.e., up to (v/c) 26 ) 
in the test-mass limit so to have sufficient analytical in- 
formation to provide a good approximation to the exact 
data. Due to recent achievements in the computation of 
the p^ix; , 0) [43], this seems a doable task. 



V. THE HORIZON FLUXES IN THE 
^-DEPENDENT CASE 

Let us now turn to discuss the ^-dependent case of 
p^ m (x; v) and in particular how to use the v = in- 
formation there. As mentioned above, Rcf. [5] argued 
that it is meaningful to hybridize the ^-dependent with 
^-independent information for the p^ n provided that the 
independence of the coefficients of their Taylor expansion 
is sufficiently mild. Following this reasoning, Ref. [5] de- 
veloped the so-called 3 +2 PN approximation to the p^: 
the ^-independent terms (up to 5PN fractional accuracy 
for the I = m = 2 multipole) added to the ^-dependent 
ones so to obtain expressions that are formally of high 
PN order, though the ^-dependence appears only in the 
low-order terms (e.g., up to 3PN for the I = m = 2 
multipole) . 

We aim here at applying the same methodology to in- 
crease the usable analytical information available for the 
pf m {x] 2/)'s. After our analysis in the test-mass limit, 
we have clearly seen that the current analytical approx- 
imation for the pf m is rather inaccurate when x in- 
creases. Luckily, we have also found that, at least for 
the £ = m = 2 case, the ^-dependence of f>22(z'i v) is 
rather mild, which suggests that an hybridization proce- 



dure can also make sense here. The unfortunate com- 
plication is that the maximum analytical information for 
the pf m is the same for both the v = and v ^ cases. 
To overcome this problem, we propose then the follow- 
ing strategy. Since the rational function representation of 
pf m (x\ 0) looks robust enough, we use it to obtain a hope- 
fully reliable estimate of the (effective) PN coefficients 
up to 4PN. We then hybridize the 1PN, ^-dependent an- 
alytical coefficient with the so-obtained test-mass 2PN, 
3PN and 4PN coefficients. Focusing on the most relevant 
I = ni = 2 multipole for definiteness, we propose to write 
/5^(£; v) as the following Taylor-expanded expression 



4- l%v + lhv 2 -2v 3 



A{l-Av + 2v 2 



„22 %,2 



„22 ,3,3 



„22 S 4 



(46) 



where the 1PN coefficient is the exact, i/-dcpcndcnt one 
coming from the analytic calculation of Eq. (40), while 
the (^-independent) c 22 coefficients with i = 2,3,4 are 
obtained from the Taylor expansion of the fitting func- 
tion to the test-mass numerical p^>(x; 0) with the 1PN- 
constrain d 22 = 1 — c 22 mentioned in the previous section. 
We finally plug Eq. (46) into Eq. (39) so to complete our 
proposal to describe the (£ = m = 2) horizon flux along 
the circular binary dynamics, including unstable orbits. 
For simplicity we consider the EOB model with the 3PN- 
accurate representation of the A(u) function discussed 
above, but there is evidently no difficulty in comput- 
ing the same quantity with the 5PN-accurate, NR-tuned 
EOB potential, i.c in the presence of the (as(^), ae(v)) 
functions of Refs. [16, 41]. 

The two panels of Fig. 4 compare and contrast the 
£ = m = 2 horizon and asymptotic fluxes, plotted versus 
x = x in the test-mass limit (y = 0, left panel) and versus 
x in the equal-mass case (y = 1/4, right-panel). The blow 
up, that for the test-mass limit happens at x hK = 1/3, oc- 
curs at x LR = 0.3625 (corresponding to tlr = 2.8456) in 
the 3PN-accurate, equal-mass case. Note that in the test- 
mass limit we present three different analytical curves: (i) 
the resummed asymptotic flux, with a 5PN-accurate ex- 
pression for p^(x; 0); (ii) the "exact" horizon flux, where 
pWix\ 0) is given by a complete fit to a rational func- 
tion, and (iii) the 4PN approximation to this latter fit 
given by Eq. (46) above. The figure indicates that this 
approximation is rather accurate up to x ~ 0.295, i.e. 
r ~ 3.4, which gives us confidence that it is meaningful 
to use it also in the equal-mass case. In addition we also 
show for completeness the i 7 ^ 1 "" 11 /^ 2 numerical points. 
In the equal-mass case the horizon flux is compared to 
the asymptotic flux computed at the 3 +2 PN approxima- 
tion from Eq. (32), notably using the f>22{x; v) of Eq. (33) 
(hybridized with 4PN and 5PN test-mass terms). 

The same procedure can be followed for the other mul- 
tipolcs, that will become more and more relevant for de- 
creasing values of v. Evidently, since the ^-dependent 
analytical information is currently limited to p22{x; v) 
only, one is obliged to use the test- mass pe m {%', 0) in PN- 
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FIG. 4. Comparison between asymptotic and horizon fluxes (divided by v 2 ) computed in the EOB approach. The left panel 
refers to the test-mass limit, while the right-panel to the equal-mass case. The left panel also highlights the difference in the 
horizon flux between the fit representation of p22(x\ 0) and its 4PN-accurate expansion. See text for further explanations. 



expanded form for the subdominant modes. Similarly, 
concerning the leading-order prefactor F^ m ' (x; v), the 
full ^-dependence is present only for the I = 2, m = 1 (see 
Eq. (14)) multipole, while for the others one can only use 
the partial independent information provided by Eq. (3). 

VI. CONCLUSIONS 

We have presented a recipe, within the EOB approach, 
to put together analytical and numerical results so to 
build an analytical expression of the horizon-absorbed 
flux of gravitational waves for circularized black-hole bi- 
naries. Our result is not limited by the weak-field, slow- 
motion approximation, but it is well-behaved also in the 
strong-field, fast-velocity regime, notably up to the last 
unstable circular orbit as defined by the EOB adiabatic 
dynamics. In this respect, our final analytical expressions 
incorporate the slow-motion results currently present in 
the literature [20, 22-24], but go beyond them adding 
new information in a resummed form. In practical terms, 
we do so by introducing a suitable factorization of the 
horizon flux, analogous to what was done some time ago 
for the asymptotic flux [5]. This factorization is such 
that the only information that needs to be computed is 
encoded in the residual amplitude correction to the hori- 
zon waveform, pf m (x; v), which plays a role analogous 
to the residual amplitude corrections to the asymptotic 
waveform, pf^ix; v) introduced in Ref. [5] and slightly 
modified here to make it well-behaved also below the 
LSO. 

In the test-mass limit [y = 0), we have computed the 
Phn( x > 0)' s numerically and we have then fitted them 
with certain rational functions so to obtain a portable 
analytical representation. The accuracy of such a fit, and 
in particular the consistency with lPN-accurate analyt- 



ical results extracted from Ref. [24] (~ 7% or smaller), 
makes us confident that it gives a reliable effective repre- 
sentation of the exact analytical pf m {x; 0) as it could be 
obtained from a high-order PN calculation able to bring 
in further corrections to the (fractional) lPN-accurate 
result of Ref. [24] . This calculation is not present in the 
literature, but we recommend it be pushed at least up to 
fractional 4PN order (i.e. 13PN in all) so to hopefully 
obtain a reasonable accuracy towards the light-ring for 
the £ = m = 2 mode. 

In the ^-dependent case, we have proposed an expres- 
sion for /5^>(x; v) given by Eq. (46) that hybridizes the 
4PN-accurate Taylor expansion of the test-mass fit with 
the lPN-accurate calculation extracted from [24]. This 
kind of hybridization can be conceptually performed also 
for the other multipoles, although ^-dependent (leading- 
order) analytical information is available only for the 
I = 2, m = 1 one, Eq. (14). The I = 3 leading-order 
expressions are implicitly present in Ref. [23] and can be 
made explicit by suitably extending the computation of 
Ref. [24]. 

The main result of this work is to provide analytical 
formulae for the horizon flux for any v that are resummed 
and that incorporate more information than the next- 
to-leading-ordcr results of Ref. [24]. The final aim of 
such expressions is to be incorporated (as a piece of the 
radiation-reaction force) in the EOB description of coa- 
lescing, nonspinning, black-hole binaries. Although for 
this kind of systems horizon absorption is expected to 
be small (since it is a fractional 4PN effect), nonetheless, 
it will introduce a correction to the phasing that might 
be relevant in ongoing comparisons with NR waveforms. 
As a preliminary exploration of this effect, we consid- 
ered the (as, a6)-dependent EOB model of Ref. [8] and 
we added a term of the form J 7 ^ = — l/(8irf2)F^ to 
the standard radiation reaction. As in Ref. [15], we se- 
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lected the values 05 = —6.37 and ag = 50 that belong 
to the region in the (05,05) plane of good agreement 
between EOB and numerical-relativity equal-mass wave- 
form found in Ref. [41] (see also Ref. [16]), and we focused 
on an equal-mass binary g = m^/ms = 1 (y = 1/4) and 
a q = 3 binary (y = 0.1875). For the g = 1 case, the 
dephasing due to our additional is ~ 5 x 10~ 3 rad 
at the light-ring crossing, and thus of the order of the 
estimated phase accuracy of state-of-the-art NR simula- 
tions, < 10~ 2 rad [28, 29]; on the contrary, for q = 3 this 
accumulated dephasing is as large as ~ 0.1 rad at light- 
ring crossing. To give more weight to this preliminary 
analysis, before embarking in the study of the spinning 
case where absorption effects are expected to be more 
relevant, it will be important to assess the accuracy of 



the fits to the test- mass pf^'s, and thus of the analyt- 
ical , in the dynamical case of BBH coalescence in 
the test-particle limit [54, 61, 62]. This work is currently 
in progress and will be presented in a separate publica- 
tion [63] 
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